Macroscopic quantum electrodynamics and duality in non-local and Onsager-violating 

media 
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We formulate macroscopic quantum electrodynamics for the most general linear, absorbing media. 
In particular, Onsager reciprocity is not assumed to hold. For media with a non-local response, the 
field quantisation is based on the conductivity tensor and the Green tensor for the electromagnetic 
field. For a local medium response, we introduce the permittivity, permeability and magnetoelectric 
susceptibilities to obtain an explicitly duality-invariant scheme. We find that duality invariance only 
holds as a continuous symmetry when non-reciprocal responses are allowed for. 
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The linear response of a macroscopic material to ap- 
plied electromagnetic fields can go beyond the scope of 
simple descriptions via electric permittivities and mag- 
netic permeabilities In particular, cross-suscepti- 
bilities naturally arise in chiral (meta-)materials Q, 
topological insulators [1] or moving media An ad- 
ditional complication arises in the latter case where the 
Onsager reciprocity [5] fails to hold. The Onsager prin- 
ciple is also violated in Tellegen media @, including 
the recently proposed perfect electromagnetic conductor 
(PEMC) that continuously interpolates between a per- 
fect conductor and an infinitely permeable material [7[. 

Chiral and anisotropic meta-materials have recently 
been discussed as candidates for repulsive Casimir forces 
[§| . Similar effects have been predicted for topological in- 
sulators [o'l and materials with Chern-Simons interaction 
(loj . Note that repulsive forces for magnetoelectric media 
do not require chiral properties; they have originally been 
discussed for dielectric plates interacting with magnetic 
ones JJj]. To implement these effects with metamateri- 
als, the anisotropic response of the medium needs to be 
taken into account [l2| . 

The impact of electric and magnetic material proper- 
ties can be studied in a systematic way by means of a du- 
ality transformation [13| . It has recently been shown that 
macroscopic QED in isotropic magnetoelectrics obeys a 
discrete duality symmetry 1J| . This has immediate con- 
sequences for dispersion forces in free space. 

In this Letter, we develop the framework for studying 
quantum optical phenomena in the presence of the most 
general linear absorbing media, including nonlocal, bian- 
isotropic and Onsager-violating materials. We show that 
for locally responding media, the electromagnetic field 
can be quantised in an explicitly duality-invariant way. 
The inclusion of Onsager-violating materials restores du- 
ality as a continuous symmetry. The successful quantisa- 
tion is in contrast to recent claims that canonical quan- 
tisation can only be performed for reciprocal media [l5| . 
Our theory lays a solid foundation for studying novel phe- 
nomena in the context of dispersion forces, Purcell effect 



and quantum light propagation. 

Field quantisation for non-local media. The constitu- 
tive relation of a linear medium can be given by Ohm's 
law in its most general form 



jin(r,i) 



dr / dVQ(r,r',T)-E(r',t-T) 



+ jN(r,0- (1) 



Here, Q(r,r',T) is the conductivity tensor and jN(r, i) 
is the random noise current required to fulfil the fluc- 
tuation-dissipation theorem pT)) as given below. Causal- 
ity requires that Q(r,r',T) = for cr < |r — r'|, in 
particular for all r < 0. Ohm's law in frequency space 
[O — doj 0{uj) + H.c] takes the simple form 

j.^(r,c.)= ydVQ(r,r',^).E(r',c^)-|-jj^(r,c^). (2) 
with 

poo 

Q(r, r', w) = 27rQ(r, r',uj)= dr e''^^Q(r, r', r) . (3) 



As a result of the causality requirement, the conductivity 
obeys the Schwarz reflection principle. 



Q*(r,r',w) = Q(r,r',-a;*). 



(4) 



Combining the constitutive relation with the Maxwell 
equations. 



p. 

V-E = =^ 



VxE - iwB = , 



V-B = 0, VxB + ^E = Aioi , 



(5) 
(6) 



one finds that the electric field obeys a generalised inho- 
mogeneous Hclmholtz equation 



VxVx 



E(r, uj) ~ ipQUj / d^r'Q(r, r', w)-E(r', w) 



i/^owjj^(i-,t^). (7) 
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Introducing the Green tensor 
VxVx ^ 



G(r,r',c.) 

-i/iow Jd^sQ{r,s,uj)-G{s,r',uj) = 5{r - r') (8) 

with G(r, r', w) ^ for |r — r'| oo, the formal solution 
to the above integro-differential equation reads 



E(r,a;) = i^o^^ / d^r'G{r,r',uj)-] {r',u!). 



(9) 



By virtue of its definition ([5]) , the Green tensor inherits 
the Schwarz reflection principle @ from the conductivity, 



G*(r,r',c.) = G(r,r', 



(10) 



As a major departure from previous treatments, we 
do not require the conductivity to obey reciprocity, i.e. 
Q^(r', r, w) = Q(r, r', w) does not necessarily hold. As a 
consequence, the Green tensor will not obey the Onsager 
principle, i.e.. 



GT(r',r,w) = G(r,r',L.) 



(11) 



will not hold in general. 

Despite this generalisation, it is still possible to de- 
rive an integral relation for the Green tensor. To that 
end, we write the Helmholtz equation ([5]) in the form 
H-G = i where (r|G|r') G(r,r',w) and (r|H|r') = 
[Vx Vx-w2/c2]^(r - r') - i^io^^Q{r,r' The Green 
operator is the right-inverse and, within any group of in- 
vertible operators, also the left-inverse of the Helmholtz 
operator, G-H = l. From this relation and its Hcrmitian 
conjugate, we find 



G-(H-Ht).Gt = Gt-G. 

In coordinate space, this relation reads 



(12) 



Mow 



J d^s J d3s'G(r,s,w)-7eeQ(s,s',w)-Gt(r',s',a;) 

= ZmG(r,r',a;). (13) 



Here, we have introduced generalised real and imaginary 
parts of a tensor field according to 

7^eT(r, r') = i [T(r, r') + T^(r', r)] , (14) 
TmT(r,r') = i[T(r,r')-Tt(r',r)]. (15) 

They reduce to ordinary real and imaginary parts for 
orthogonal tensor fields with T^(r',r) = T(r, r'). The 
integral relation generalises the result from Ref. [ll| 
to the case where Onsager reciprocity does not hold. 

Returning to the electromagnetic field, an explicit 
quantisation is achieved by specifying the commutator 



T(r,w),jj^(r',a;') 



TZeQir, r\uj)6{uj -uj'). (16) 



The fact that the right-hand side is a Hermitian tensor 
field guarantees the consistency of this commutation rela- 
tion. Introducing the ground state |{0}) of the medium- 
field system according to jj^(r,a;)|{0}) = 0, the currents 
satisfy the fluctuation-dissipation theorem 

({AjN(r,w),Ajj,(r',a.')}) 

= - Xm[icjQ(r, r', uj)]S{uj - uj') . (17) 

TT 

Combining Eqs. ([9]) and (|16p . one flnds that the fluctu- 
ations of the electric field are also consistent with the 
fiuctuation-dissipation theorem, as required: 

({AE(r,w),AEt(r',c.')}) 

^-Im\^ouj^G(r,r',u})]S(u}~u}'). (18) 

To verify the canonical equal-time commutation rela- 
tions, we introduce the vector potential for the electro- 
magnetic field in the Coulomb gauge, A = E^/(iaj) (_L: 
transverse part). Using Eqs. ([9]) and ([T6|) . one finds 



[E(r,c.),At(r',^')] 



ZmG-L(r,r',a;)(5(w - w') (19) 



and hence 
[E(r),A(r')] 

= ^ y'°°d....[G^(r,r',a.) + ^GT(r',r,a.)] (20) 

where the Schwarz refiection principle pUj) has been 
used. Closing the integration contour in the upper 
half of the complex uj plane and using the asymptote 
(a;^/c^)G(r, r', w) — —6{r — r') for w — ?• oo, one finds the 
canonical commutation relation from free-space QED, 



[E(r),A(r')] =-5^(r-r') 



(21) 



Next, we introduce bosonic creation and annihilation 
operators f^, f according to the prescription 



jN(r,cj) = 



dVR(r,r',w) •f(r',w) 



(22) 



where R is a square root of the positive definite tensor 
field 7^eQ, 



J d^sR(r,s,w)-R'l'(r',s,w) = 7leQ(r, r',uj). 



(23) 



Together with Eq. (|T6l) . this ensures bosonic commuta- 
tion relations, [f(r,a;),f1'(r',w')]='5(r-r')<5(a;-w')- The 
Hamiltonian of the body-field system is then 

HF=Jd^rJ^ dLunuj{^r,uj)-i{r,uj). (24) 

It leads to i{r,uj,t) = f (r, aj)e^'"*, hence Maxwell's 
equations for the electromagnetic-field operators in the 
Heisenberg picture are valid by construction. 
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Field quantisation for local bianisotropic media. For 
media with a spatially local response, it is convenient to 
cast the inhomogeneous Maxwell equations dH) into the 
alternative forms 



V-D = 0, VxH + iwD = 



(25) 



with 



D = £oE + P, H=— B-M. (26) 
Mo 

The polarisation and magnetisation fields respond lin- 
early to the electric and magnetic fields {Zq = /io/eo ), 

P - eo(£-|-/x-i-C-l)-E + Zo-i|.^-i.B + PN, 

(27) 

M = VV"'-C-E + Mo'(l - At"')-B + MN, (28) 

where e(r,a;) is the medium's permittivity; fi{r,uj) its 
permeability; ^(r,w) and C(r,'^) its magnetoelectric sus- 
ceptibilities; P^ and Mjj denote the noise polarisation 
and magnetisation. Combining Eqs. (|26 |) -([28 l) . the con- 
stitutive relations can be given in the more familiar form 

D = £o£-E + c-i^-H + Pn + c-i^-Mn, (29) 

B = c"^C-E + moM-H + ^oM-Mn- (30) 

Combining Eqs. ([25]), ([291) and dSO]), we see that they 
are a special case of Eq. ([7]) with 

Q(r,r',a;) = (i^o'^)-i V x (/x^i - 1) •<5(r-r') x V' 

+ Zq-^V X /i-i-C-<5(r-r') + Z^^$-fi-^-5{r~r') x V' 
-ieoio{e-^-H-^-C-\)-6{r~r'), (31) 

jj^ = -ic^EN + VxMN- (32) 
The Green tensor for the electric field dHl) solves 



Vx/x^^-Vx -— Vx/x^^-C- +— |-/x"^-Vx 

c c 



^(e-^/x-^-O 



G(r,r',a;) = 5(r-r'). (33) 



The commutation relations for Pn and Mn can be 
deduced by substituting Eqs. dSJ) and into Eq. ([TC)). 



[PN(r,^),PMr-,w' 



Xni(£-^-/x-i-C)(5(r - r')<5(w - w'), (34) 
[PN(r,^),MUr',^')] 

h 



(C^-M-it-^.;,-i)5(r-r')<5(c.-^'), (35) 



27riZo 
[MN(r,^),MUr',c.')] 

= —ImiJi^H{v - v')5{u - Lo'). 



(36) 



We introduce bosonic creation and annihilation operators 
(ix (r, w) , il (r', u')] = 5xyS{v - v')5{u - c^') (A, A' = e, m) 
according to 



■PN(r,c.)\ /ft^Yfe(r,c.) 

,MN(r,w)y V TT Vfm(l-,'^) 

where the (6 x 6)-matrix 7?. is a root of 



(37) 



-It.tt 



2iZn 



2iZo_ 
Im/x ^ 

Mo 



(38) 



The Hamiltonian of the body-field system is again 
quadratic and diagonal in the bosonic variables, 



E 



dLuf^{l{r,Lj)-h{r,Lo). (39) 



Note that Eqs. imply a separation of the in- 

ternal current density into electric and magnetic parts, 
i = — iwP-f VxM. This separation and the resulting 
explicit field quantisation is not unique [17| . 

Duality invariance. Introducing dual-pair notation 
(E"r,ZoH"r)"^, (ZoD"^,B"r)"r, we may write the MaxweU 
equations ([5]) and ([25)) in the compact form 



B 

1 

-1 






ZoR 

B 



(40) 
(41) 



The constitutive relations (12^]) and ([50]) read 



ZoD 
B 



E 

ZoH 



1 € 
n 



^oPn 
MoMn 



(42) 

Maxwell's equations are invariant under duality trans- 
formations 



® 



cosf 
~ sin^ 



smt 
cos( 



(43) 

because the symplectic matrix in Eq. ()4ip commutes 
with D{9). From the constitutive relations, we find 
transformed medium response functions (£,^,^,/x)^® = 
Vie)ie,tC,f^V with 



v{e) = 

( COs2 d 

~ sin B cos i 



\ sm 

as well as 

1 I 
/2 



sm f cos f 
2 a 



sm V cos ( 



— sm V 
sin 9 cos ( 



sm W cos W 
- sin^ e 
cos^ 
— sin 9 cos ( 



sm y 
sin 6 cos 
sin cos 

cos2 6' 



■^oEn 
MoMn 



= D{e) 



1 





^oE] 



■N 



MoMn 



(44) 
\ 

1 

/ 

)■ 

(45) 
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It is worth discussing a few special cases of bian- 
isotropic media, their characteristic features and be- 
haviour under duality transformations: 

• Isotropic media (£ = el, fi^fil, ^ = ^ = 0): Onsager 
reciprocity ([TT|) holds; Pn and M]^ commute; gen- 
eralised real and imaginary parts reduce to ordinary 
ones; discrete duality symmetry. 

• Biisotropic media {e — = ^ = C ~ CO- 
generalised real and imaginary parts in Eqs. (j34p 
and p6|) reduce to ordinary ones; continuous dual- 
ity symmetry. 

• Anisotropic media = C = Pn and mJj com- 
mute; discrete duality symmetry. 

• Reciprocal media (e^ — e, = — /^^ = A*): 
Eq. (jlip holds; generalised real and imaginary parts 
reduce to ordinary ones; discrete duality symmetry. 

Here, discrete duality symmetry means that the rotation 
angle is restricted to values 9 = mT/2 with nGZ. 

To distinguish reciprocal magnetoelectric susceptibili- 
ties from non-reciprocal ones, one commonly writes ^ — 
— i^^ and C = X + i**- Here, the chirality tensor 
K ~ (C"-^^) / (2i) represents the reciprocal magnetoelectric 
response; and the non-reciprocal magnetoelectric tensor 
X = (C+€^)/2 vanishes for a reciprocal medium. 

To derive transformation laws for the Green tensor, we 
combine Eqs. ©, Q, ([Ml), and ^ to write 

(A) 

(47) 

where integration over the second position variable is 
implied and we have introduced the shorthand nota- 
tions Gee = {iuj/c)G{iuj/c), Gem = {iLO / c)G X V', G^e = 

V X G(ia;/c), and G,„m = V x G x V'. The trans- 
formed tensors follow upon applying duality transfor- 
mations on both sides of the equation. In general, 
the result is very complex due to the presence of the 
medium response functions in Eqs. psj) and ([Tf]). In 
the special case of r and r' being in free space, one ob- 
tains the simple result (Gee, Gem, G^e, Gmm + ^)^® = 
I?(6')(Gee, Gem, Gme, Gmm + ^)^- The Green tensors 
then transform like the medium response functions with 

Eq. m. 

Conclusions. Starting from the general Ohm's law, 
we have quantised the electromagnetic field in the pres- 
ence of non-local, non-reciprocal media which satisfies 
(i) the canonical commutation relations from free-space 
QED; (ii) the linear fluctuation-dissipation theorem; and 
(iii) the macroscopic Maxwell equations. Key feature 



of the scheme is the symmetrisation of tensor fields via 
generalised real and imaginary parts, which is necessary 
whenever Onsager reciprocity does not hold. 

For a local bianisotropic medium, we have shown that 
quantisation can alternatively performed by introducing 
permittivity, permeability and magnetoelectric suscepti- 
bilities. When the latter do not vanish, the noise polari- 
sation and magnetisation do not commute. We have ex- 
plicitly determined the behaviour of the fields, response 
functions and Green tensor under duality transforma- 
tions. The full continuous transformation group applies 
for bianisotropic and biisotropic media, but reduces to 
a discrete symmetry for anisotropic, isotropic and/or re- 
ciprocal media. 

The scheme lays the foundation for exact studies of 
quantum phenomena such as dispersion forces, (Forster) 
energy transfer or environment-assisted molecular tran- 
sition rates in the presence of motion or novel media with 
chiral or non-reciprocal properties. For example, it will 
facilitate the design of materials that show Casimir repul- 
sion, and guide towards new schemes to test CP-violation 
in atoms. 
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Physical Sciences Research Council (EPSRC). We would 
like to thank F. H. Hehl for discussions. 



* Electronic address: |s.buhmann@imperial.ac.uk| 
[1] L.V. Lindell, A.H. Sihvola, S.A. Tretyakov, and A.J. Vi- 
itanen. Electromagnetic Waves in Chiral and Bi-Isotropic 
Media (Artech House, Norwood, MA, 1994); F.H. Hehl 
and Y.N. Obukhov, Foundations of classical electrody- 
namics (Birkhauser, Boston, 2003); Phys. Lett. A 334, 
249 (2005). 

[2] A. Ishimaru, S.W. Lee, Y. Kuga, and V. Jandhyala, IEEE 
Trans. Antennas Prop. 51, 2550 (2003); V. Yannopapas, 
J. Phys.: Condens. Matter 18, 6883 (2006); M. Thiel, 
M.S. Rill, G. von Freymann, and M. Wegener, Adv. Mat. 
21, 4680 (2009). 

[3] L. Fu, C.L. Kane, and E.J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[4] J.F. McKenzie, Proc. Phys. Soc. 91, 532 (1967). 

[5] L. Onsager, Phys. Rev. 37, 405 (1931). 

[6] B.D.H. Tellegen, Philips Res. Rep. 3, 81 (1948). 

[7] I.V. Lindell and A.H. Sihvola, J. Electromagn. Waves and 
Appl. 19, 861 (2005). 

[8] R. Zhao, J. Zhou, T. Koschny, E.N. Economou, and 
CM. Soukoulis, Phys. Rev. Lett. 103, 103602 (2009). 

[9] A.G. Grushin and A. Cortijo, Phys. Rev. Lett. 106, 
020403 (2011). 

[10] V.N. Marachevsky and Y.M. Pis'mak, Phys. Rev. D 81, 
065005 (2010). 

[11] C. Henkel and K. Joulain, Europhys. Lett. 72, 929 
(2005); M.S. Tomas, Phys. Lett. A 342, 381 (2005). 

[12] F.S.S. Rosa, D.A.R. Dalvit, and P.W. Milonni, Phys. 
Rev. Lett. 100, 183602 (2008). 

[13] L. Silberstein, Ann. d. Phys. 327, 579 (1907); G.Y. 
Rainich, Trans. Am. Math. Soc. 27, 106 (1925); C.W. 



5 



Misner and J.A. Whcclcr, Ann. Phys. 2, 525 (1957). 
[14] S.Y. Buhniann and S. Schccl, Phys. Rov. Lett. 102, 

140404 (2009). 
[15] T. Philbin, New J. Phys. 13, 063026 (2010). 
[16] C. Raabe, S. Scheel and D.-G. Welsch, Phys. Rev. A 75, 



053813 (2007). 

[17] D.B. Melrose, Quantum Plasmadynamics (Springer, New 
York, 2008). 



